Single magnetic molecule between conducting leads: Effect of mechanical rotations 
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We study spin-rotation effects in a magnetic molecule bridged between two conducting leads. Dy- 
namics of the total angular momentum couples spin tunneling to the mechanical rotations. Landau- 
Zener spin transition produced by the time-dependent magnetic field generates a unique pattern of 
mechanical oscillations that can be detected by measuring the electronic tunneling current through 
the molecule. 

PACS numbers: 75.50.Xx,85.65.-|-h,75.45.-|-j,85.75.-d 
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In the last years a significant experimental effort has 
been made to measure electronic current through a sin- 
gle magnetic molecule bridged between conducting leads. 
These studies have been driven by possible applications of 
single-molecule magnets in spintronics, as well as by the 
hope to use magnetic molecules as qubits [I] . Heersche et 
al. observed striking voltage dependence of the current 
through a single- molecule magnet Mni2, with a complete 
current suppression and excitations of negative differen- 
tial conductance on the energy scale of the anisotropy 
barrier [3]. Jo et al. measured the magnetic field depen- 
dence of the electron tunneling spectrum in a transistor 
incorporating a Mni2 molecule [3]. Henderson et al. ob- 
served a Coulomb blockade effect by measuring conduc- 
tance through a single-molecule magnet Mni2 Voss 
et al. conducted experiments and developed theoretical 
model for the dependence of the tunneling current on the 
orientation of the Mni2 molecule 0j. 

Various aspects of the electronic transport through 
magnetic molecules have been investigated theoretically. 
The effect of the exchange coupling between spins of con- 
ducting electrons and the spin of the Mni2 molecule has 
been studied by G.-H. Kim and T.-S. Kim [6;']. Elste and 
Timm [Tj developed a model for the Coulomb blockade 
in a transport through a single magnetic molecule weakly 
coupled to magnetic and nonmagnetic leads. They also 
studied the possibility of writing, storing, and reading 
spin information in memory devices based upon single- 
molecule magnets Q . Kondo effect in transport through 
a single-molecule magnet strongly coupled to metallic 
electrodes has been investigated by Romeike et al Q. 
The effect of spin Berry phase on electron tunneling has 
been studied by Gonzalez and Leuenberger ■ Misiorny 
et al. investigated magnetic switching of the molecular 
spin by a spin-polarized current [Il | , as well as tunneling 
magnetoresistance [12|. Cornaglia et al. have studied 
the effect on the transport of a soft vibrating mode of 
the molecule 13|. First principle DFT calculations of 



the electron transport through a Mni2 single-molecule 
magnet and of spin-filtering effect have been performed 
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by Barraza-Lopez et al. [IJ]. There also has been some 



effort to compute Josephson current through a magnetic 



molecule coupled to superconducting leads [15| . 

Recently two of the authors have demonstrated 
that quantum states of a magnetic molecule that is free 
to rotate are different from quantum states of a magnetic 
molecule whose position is fixed in a crystal. This ef- 
fect arises from the conservation of the total angular mo- 
mentum (spin -I- orbital). Somewhat intermediate situa- 
tion occurs when a molecule is bridged between two elec- 
trodes. While such a molecule is not free to rotate, its ori- 
entation can still vary due to the spin-rotation coupling. 
In this Letter we study the Landau-Zener transition be- 
tween quantum states of the molecule caused by a field 
sweep. We show that coupled dynamics of the molecular 
spin and the mechanical rotation of the molecule leads 
to the modification of the Landau-Zener dynamics. It 
also results in a specific pattern of mechanical oscilla- 
tions of the molecule. Since electron tunnelling through 
the molecule should strongly depend on its orientation, 
this effect should be seen in the tunneling current. 

We consider geometry depicted in Fig. [TJ The mag- 
netic anisotropy axis of the molecule, Z, is perpendic- 
ular to the direction of the transport current Y. The 
mechanical motion of the molecule is restricted by its 
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coupling to the leads. For simplicity we consider only 
small torsional oscillations about the Z-axis that we de- 
scribe by the angle of rotation (j>. This is justified by the 
fact that quantum tunneling of the molecular spin be- 
tween two opposite directions along the Z-axis changes 
only the Z-component of the angular momentum, thus 
generating the ^-component of the torque. In our ap- 
proximation, small rotation of the molecule due to this 
torque produces negligible deformation of a rigid mag- 
netic core of the molecule. On the contrary, the tunnel- 
ing current through the molecule must have exponential 
dependence on (j). We assume that electrons in the leads 
have only marginal effect on the quantum spin states of 
the molecule. 

The main part of the Hamiltonian of the molecule has 
three terms: 



Hz + Hs 



Here 



(1) 



(2) 



describes mechanical rotations about the Z-axis. Iz is 
the corresponding moment of inertia and u>r is the fre- 
quency of free torsional oscillations of the molecule due 
its coupling to the leads. Operator of the mechanical 
angular momentum, = —id/d(j), satisfies the commu- 
tation relation 



, L.] = i. 



(3) 



Interaction of the molecular spin S with the external 
magnetic field, B applied along the Z-axis, is described 
by the Zeeman term in Eq. ([T]): 



Hz = -g^J.BSzBz 



(4) 



with g being the gyromagnetic factor and fiB being the 
Bohr magneton. Operator Hs in Eq. ([T|) is given by 

Ulii 



Hg — RHaR 



(5) 



where Ha is the crystal-field (magnetic anisotropy) spin 
Hamiltonian and 



R 



is the operator of rotation in the spin space. 
The general form of Ha is 

Ha = H\\+H^, 



(6) 



(7) 



where commutes with Sz and H± is a perturbation 
that does not commute with Sz- The existence of the 
magnetic anisotropy axis Z means that the | ± S) eigen- 
states of Sz are degenerate ground states of H\\. Op- 
erator H± slightly perturbs the | ± S) states, adding 



to them small contributions of other \ms) states. We 
shall call these degenerate normalized perturbed states 
\ip±s)- Physically they describe the magnetic moment of 
the molecule looking in one of the two directions along the 
anisotropy axis. Full perturbation theory with account 
of the degeneracy of Ha provides quantum tunneling be- 
tween the \ip±s) states. The ground state and the first 
excited state become 



They satisfy 



(8) 



(9) 



with E^—E^ = A being the tunnel splitting. Expressing 
iV'is) via according to Eq. ([5]), it is easy to see from 
Eq. 11]) that 

{'^±s\HA\i^±s) = 0, (V'-sI^aI^s) = -A/2 . (10) 
This gives 

A 



(V'±5|^5|^±5) = 0, {ll}^s\Hs\'lp±s) 



^±2iS<t> 



(11) 

for the matrix elements of Hg. 

Since the low-energy spin states of the molecule are 
superpositions of \i)±s)^ it is convenient to describe such 
a two-state system by a pseudospin 1/2. Components of 
the corresponding Pauli operator cr are 



fJx = \i'-s){tps\ + \^s){'>P-s\ 

(jy = i\%l)_s){'^s\- i\i's){'4>-s\ 

(^z = |V's)(V's| - |'/'-s)(V'-s| ■ 



(12) 



The projection of any operator H onto \i)±s) states is 
given by 

^ (l^™|i?|V'n)IV'm)(V'n|- (13) 
m,n— ±5 

Noticing that 

Sz\^±s)'^Sz\±S)^±S\^±s). (14) 

it is easy to project Hamiltonian (IT|) onto t\}±s- Simple 
calculation yields 

^cff = ^rot - ^VKa^ - ^A [cos {2S4>) ^ sin (2S^) Uy] 

(15) 



(16) 



where 

Heff = Wgz + A cos {2S(j)) + A sin {2S(j)) Gy (17) 



3 



and 



(18) 



Note that for a non-rotating magnetic molecule the 
effective Hamiltonian has the form 



Its eigenvalues are 



1. 



1 



-AcTt 



with 



(19) 



(20) 



(21) 



Using the Heisenberg equation of motion for an oper- 
ator A{t), 



Ait) = - [H,s,A{t) 



(22) 



one obtains 



nt^ = -IzOJrCl) - SA [sin {2S(p) cr^ - cos {2S(j)) ay] 



—r- , a =-[a X Hoff 
L n 



(23) 



Elimination of provides the following system of equa- 
tions: 



position, = 0. In the absence of such a torque the to- 
tal angular momentum of the molecule alone would be 
conserved. 

Equations l|24p - (P7)l are operator equations. To 
obtain numerically tractable equations, these operator 
equations should be averaged over quantum states of the 
system. If one decouples the quantum averages as 



(sin (ip) 



(31) 



in the spirit of the mean-field approximation, one obtains 
classical-like equations of the same structure as above. 
Without such a decoupling the equations for the Heisen- 
berg operators are useless and one has to go back to the 
Schrodinger equation for the whole system consisting of 
coupled spin and mechanical subsystems. While the de- 
coupling cannot be justified in the general case, we notice 
that the cross terms containing cr and (f) become small 
when (if) is small. Since ip = 2S(j) this condition for 
large S is stronger than the condition of small oscilla- 
tions, <§i 1. Nevertheless, even for S = 10, which is 
the case of a Mni2 molecule, it is likely that rotations of 
the molecule in the geometry depicted in Fig. [1] will still 
satisfy the condition \(p\ ^ 1. In this case Eqs. ([M]) - 
(P7)) can be treated as classical. 

The behavior of the system depends on the dimension- 
less magneto-mechanical constant of a free molecule. 



5 = 



hujc 
"a" 



2q2 



2h^S 



(32) 



and 



Here 



and 



h&x = Way — A sin (ip) az 
h&y — A cos (ip) az — Wax 
h&z — A sin (p) ax — A cos (p) ay 



p + (JIp = Lul [cos {(p) ay — sin [p) ax 



p = 2S(j) 



y/2S^A/h. 



(24) 
(25) 
(26) 



(27) 



(28) 



(29) 



Combining the above equations it is easy to see that 
they satisfy 



— (?),L2 + nSaz) = -IzUJr<i 



(30) 



The left-hand side of this equation is the time deriva- 
tive of the Z-component of the total angular momentum, 
h{Lz + Sz), while the right-hand side is the elastic torque 
due to the coupling to the leads (see Fig. [T|) . This torque 
acts such as to return the molecule to its equilibrium 



and the resonance parameter, 

huJr 



(33) 



For a linear field sweep, W = vt, another relevant dimen- 
sionless parameter is 



2hv ' 



It determines probability. 



P = e- 



(34) 



(35) 



of staying in the initial ■f/'-s state in the standard Landau- 
Zener problem. 

Using the above parameters and dimensionless vari- 
ables 



Vt hvT TTT 



"=T ^=A=A^=27' 



(36) 



Eqs. ([M]) - ((?7|) for \p\ ^ 1 can be re- written in the form 



a^ = way - paz 



Oz = y^'^x ^ <^y 



(37) 
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and 



(38) 



where prime means d/dr. Notice that in the last equation 
we introduced a term jip' that describes damping of the 
mechanical oscillations of the molecule. 




FIG. 2; Typical time dependence of the rotation angle (p. The 
inset shows the fine structure of the oscillations. 



The above equations have been solved numerically us- 
ing Mathematica. With the Mni2 molecule in mind, 
we compute the time dependence of the angle of twist, 
4> = ip/l2S), for = 10, ^ = 10^ and r = 10^ (see dis- 
cussion below). For the illustration purpose, the value of 
the damping constant 7 = 0.1 has been chosen to provide 
the finite duration of the oscillations. Fig. [5] shows the 
results for a relatively fast field sweep corresponding to 
e — 0.1. For a slow sweep (e 3> 1) the amplitude of the 
oscillations is significantly smaller. 

Simple arguments allow one to understand the behav- 
ior shown in Fig. [21 Oscillations of </> are excited when the 
time-dependent distance between the spin levels given 
by Eq. (|2ip , coincides with the frequency tOr of the me- 
chanical oscillations of the molecule. For large er this 
happens at r = 2er/7r in excellent agreement with Fig. 
[2j Oscillations continue at the frequency LOr until they 
are completely damped due to the finite 7. To obtain 
the dependence of the initial amplitude of oscillations on 
the parameters the approximate solution of Eqs_^ ( 1371 ) for 
ay{T) at small ip, e < 1 and r > 1 can be used 



ay{T) = -2^6-^(1-6-^) 



sm 



4e 



This gives an approximate solution of Eq. 



(39) 



at 7 = 0: 



-Im 



dr'e- 



(40) 



The result of the integration can be expressed in terms 
of trigonometric and error functions. While the formula 



is rather cumbersome it provides a simple dependence of 
the final amplitude of the undamped oscillations on the 
parameters: 



2ns 

IzUJr 



Vee-<^ (1 - e-^) . 



(41) 



It reaches maximum at e — 1.45. This expression is in 
excellent agreement with numerical results. It can be 
used as long as 2S4>max ^ 1- 

Most of the existing experiments on transport through 
individual magnetic molecules have been done with Mni2 
acetate. The moment of inertia of the Mni2 molecule 
is in the ballpark of 10~^''g-cm^. The tunnel splitting 
A depends strongly on the transverse anisotropy, H± in 
Eq. I?]). For the molecule bridged between two leads it 
may differ drastically from that in a Mni2 crystal. With 
Iz ^ lO^'^^g-cm^ our choice of (5 ~ 10^ in Fig. ^ cor- 
responds to A/h ~ lO^s"^. The resonance frequency of 
the rotational oscillations of the molecule must be much 
higher. The choice of r = 10^ in Fig. ^ corresponds to 
u!r — lO^s^^. The mechanical oscillations of the molecule 
begin at t = 2eh'^u)r/{'^A'^) after the spin states -0-5 and 
ips cross due to the field sweep. For the choice of pa- 
rameters used in Fig. ^ this time is of order lO^'^s. 



Spin oscillations given by Eq. (|39p must dissipate on a 
longer time scale for the mechanical oscillations to be ob- 
servable. In a Mni2 crystal, dissipation is dominated by 
phonons. For a setup shown in Fig. [1] phonon processes 
must be suppressed, making spin relaxation times of or- 
der lO^^'s quite reasonable. For A/h ^ lO^s^^, the fast 
field sweep, e < 1, corresponds to a few kOe per second, 
which is also reasonable. Magnetic molecules other than 
Mni2 acetate may prove to be even better candidates for 
such an experiment. 

If a tunneling current flows through the setup depicted 
in Fig.[Tl the amplitude of the current should be affected 
by the mechanical oscillations of the molecule. Same as 
in tunneling microscopy, this effect should be detectable 
due to high sensitivity of the electron tunneling rate to 
the orientation of the molecule [s*]. Following the field 
sweep, the current should acquire an oscillating compo- 
nent similar to that shown in Fig.[2l Such an experiment 
is practicable and it would be of great interest as it would 
prob quantum dynamics of spin in individual magnetic 
molecules. 

This work has been supported by the U.S. National 
Science Foundation through Grant No. DMR-0703639. 
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